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We derive an Abelian-like Ward identity in color superconducting phase and calculate vertex 
corrections to the color superconducting gap. Making use of the Ward identity, we show that 
subleading order contributions to the gap from vertices are absent for gapped excitations. 



. I. INTRODUCTION 

O 

Quark matter at large quark chemical potential (i is a weekly coupling system because momenta exchanged in the 
interaction between quarks near the Fermi surface is of order //, which makes the coupling constant g small due to the 
+^ ■ property of asymptotic freedom in quantum chromodynamics (QCD). In this case, the dominant interaction between 

Otwo quarks is one-gluon exchange, which is attractive in the color-antitriplet channel. Consequently, at sufficiently 
low temperatures, the quark Fermi surface is unstable with respect to the formation of Cooper pairs [3, [H which 
leads to the so-called color superconducting (CSC) state 0, i, I, i, 0, 0, i, EE El El, (for reviews, see, e.g. 

i, m e m ee m, m, m m ) ■ 

In a superconductor, exciting particle- hole pairs costs at least an energy amount 2</>o, where <f>o is the value of the 
' (— i ' | superconductor gap parameter at the Fermi surface for T — and can be computed from a gap equation derived under 
O ■ mean-field approximation which involves one-gluon exchan ge a nd bare quark-quark-gluon vertex. Schematically this 
I | gap equation can be written in the form [IE [22, HE HE 125L l26| 
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For small value of the QCD coupling constant, g <C I, the solution is 

c 



(1) 
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°° . , . , 

VO ■ 0o-26 M exp — )[l + 0{g)} (2) 

The first term in Eq. JTJ) contains two powers of the logarithm \a.{ji / 4>o) , one is the same as in BCS theory fl|, 2f and 
the other from the exchange of almost static magnetic gluons, which is a long-range interaction 0, EH, El, \m Thc 
weak coupling solution (|2|) implies that this term contributes to the gap equation at the order O(l). We call this 
term the leading order term. The value of the coefficient £ determines the constant c in Eq. ([2]) . The second term 
in Eq. ([T]) contains subleading contributions of order O(g) to the gap equation, characterized by a single power of 
the logarithm ln[ii/</>o) ~ 1/g. Part of it arises from the exchange of non-static magnetic and static electric gluons 
0, [S3, [H, [H, [25| . Another source is the quark self-energy correction [23] . The coefficient j3 in Eq. fT} determines 
the constant b in Eq. The term is called subleading one. The third term in Eq. |T]) summarizes sub-subleading 
contributions of order 0(g 2 ) with neither a collinear nor a BCS logarithm. It was argued in Ref. [H EE H3] that 
at this order gauge-dependent terms enter the QCD mean field gap equation. In Coulomb gauge the authors of Ref. 
[28| showed that the gauge-dependent contribution arising from the gluon propagator appear at this order when the 
momentum arguments of the gap are put on the quasi-particle mass shell. In covariant gauge one can see that the 
gauge dependence shows up at subleading order and brings an additional factor exp(3/2£) [llj to prefactor b. However, 
the gap parameter is in principle an observable quantity on the quasi-particle mass-shell, and thus gauge-independent. 
Therefore, the naive mean-field approach to the gap equation with bare qqg vertices is not enough to guarantee the 
gauge independence even at the subleading level. 

A general way to study gauge independence is to make use of Ward identities. This approach has been frequently 
used to show the gauge independence of physical collective excitations in thermal gauge theories, like hot QCD 
[2^. [30L [3ll |32| . However, we need to use Nambu-Gorkov (NG) basis in color superconducting phase, therefore it is 
desirable to derive a Ward identity in NG basis with diquark condensates. Recently Gerhold and Rebhan (33j used 
generalized Nielsen identities to give a formal proof that the fcrmionic quasiparticle dispersion relation in a color 
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superconductor are gauge independent under the assumption that the 1PI part of variation induced by that of gauge 
fixing function in the effective action has no singularities coinciding with those of quark propagator. We have provided 
another proof of gauge independence of 2SC gap by deriving a generalized Ward identity from QCD with diquark 
condensate and by applying it to gap equation [34j . 

In this paper, we present an investigation of the vertex contributions in the gap equation. The calculation is done in 
Nambu-Gorkov (NG) formalism in super-phase with diquark condensates. We will show that this method is equivalent 
and a good alternative to that used in Ref. [l2j based on four-fermion scattering amplitudes. We found that there is a 
similar cancellation as in normal phase between the Abelian and triple-gluon vertices, which leads to an Abelian-likc 
Ward identity in the NG form except that an additional term appears in the super-phase. With this Ward identity, we 
finally show that the contributions from vertices to the gap equation are free of subleading terms for gapped modes. 

In this paper four- momenta are denoted by capital letters, = (fc , k), with k being a three-momentum of modulus 
|k| = k and direction k = k/fc. For the summation over Lorentz indices, we use a notation familair from Minkowski 
space, with metric g^ v = diag(+, — , — , — ). For simplicity and without ambiguity we always write Lorentz indicies as 
subscripts. 



II. SETTINGS 



In this section, we will give some preparation knowledge and conventions necessary to the calculation. Since we 
are concerned with the super-phase, it is quite natural to work in NG basis. We will see that it is very convenient to 
describe the mass-shell condition for quasi-particles in NG basis. In this paper we choose a special case for convenience, 
the color superconducting phase with two flavors (2SC). The calculation can be extended to other phases. We work 
in zero temperature. 

In NG basis the quark fields read 

* = (| c ), * = &>) (3) 

where conjugate fields are given by tp c = Cip and ip c = ip T C with charge conjugate matrix C = ij 2 jo- The quark 
propagator inverse is 



E21 E 22 J ' W 

where we use boldface letters to denote NG matrices. The free part is given by 5 n (K) = If + fijo and S'22 C^O = 
H ~ ^70- I n super-phase, the off-diagonal elements of selfenergy S are proportional to the diquark condensate or the 
gap parameter, 

E21 = J 3 T 2 75</> 6 A e 

E i2 = -J 3 T 27 5<f *A~ e , 

where the color part is chosen as {Js)ij = ici3j, the third generator of SO(3), to incorporate the pairing in the 
anti-symmetric channel between the red and green fundamental colors. The flavor part T2 is the second Pauli matrix 
representing the anti-symmetry in flavor space. The appearance of 75 implies that we only consider even-parity 
pairings. Then the quark propagators are given by finding the inverse of S _1 (AT) in Eq. ([¥]), 
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where the repetition of indices implies summation if not indicated explicitly. Note that there are two branches 
of excitations, one is gapped denoted by i = 1, the other is gapless denoted by i = 0. Then we have 4>i = 4> 
and </>o = 0. The quasi-particle energies are e qjje = y/ (q — efi) 2 + \<ftf\ 2 . The color projectors are L\ = Jf and 
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Ln = 1 — Jo, c orresponding; to the gapped and gapless mode respectively. The quark wave- function renormalization 
[23l [35l . [36l. l37l [H, [H, H3| constant Z(q ) results from the diagonal part of the quark selfenergy, En and E 2 2> an d 
is given by Z(q ) = 1 — (g 2 / 18n 2 ) \n(M 2 / q^) with M 2 = Nfg 2 fi 2 /(2n 2 ), a scale characterizing Debye or Meissner 
screening. For two loop corrections to the gap equation in NG basis, Z(qo) can be neglected because its contribution 
is beyond sub-subleading order. The bare quark-gluon vertex is 



i(0)o _ 



T a q 

-T al 



1„ = T a 7 M 



(6) 



Here we write it in a special way with color and Dirac part separated. 

In this paper, we choose covariant gauge, the hard dense loop (HDL) propagator for gluons reads 



pL 



p 2 + ill p 2 + n t * p 4 

where n; and n t are HDL self-energies of longitudinal and transverse gluons respectively and are given by 



(7) 
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The projectors P^ v and PT V are 



Piai — 9lJ.v 
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In the Landau damping region, |po| p, n/(P) w P 2 M 2 /p 2 and n t (P) w — (7r/4)M 2 |w p |/p, where we used lo p = ipo 5 
the energy in Euclidean space or the Mustubara frequency in finite temperature. Here the subscript of tu p only means 
it belongs to the four-momentum P = (pa, p) = (— iu> p , p). 

In NG basis, the mass shell condition for quasi-particles can be expressed as 

S" 1 ^)*^) = 0, 

*(if „)S- 1 (if on ) = 0, (10) 

where K on — (e£ if k) denote on-shell momenta and ^(K on ) are on-shell wave functions. Note that these equations 
are in matrix form. Hereafter we suppress the subscript of K on for simplicity of notations if without ambiguity. 



III. CORNWALL-JACKIW-TOMBOULIS FORMALISM 



The gap parameter in superconducting systems is not accessible by means of perturbation theory: one has to 
apply non-perturbative, self-consistent, many-body resummation techniques to calculate it. For this purpose, it is 
convenient to employ the Cornwall- Jackiw-Tomboulis (CJT) formalism [41(. The first step is to add to the action 
local source terms coupled to each field and the bi-local source terms coupled to each pair of same fields. One then 
performs a Legendre transformation with respect to all sources and arrives at the CJT effective action [4l|. The 
expectation values for the one- and two-point functions of the theory are determined from the stationarity conditions 
which give Dirac equations for quarks, Yang-Mills equations for gluons and Dyson-Schwinger equations for gluon and 
quark propagators. 

In order to include all subleading contributions to the gap equation arising from the off-diagonal components of 
Dyson-Schwinger equations for quarks, we should consider all 2PI diagrams up to two-loops. There are three 2PI 
vacuum graphs as shown in Fig. (fTJ). Note that all these diagrams are in NG basis. Cutting one quark leg leads 
to all one-particle-irreducible (1PI) diagrams contributing to the quark selfenergy in Dyson-Schwinger equations, see 
Fig. {TJ(al,bl,cl). As we are mainly concerned with the gap equation, we expand the 21-components of all these 1PI 
diagrams in terms of their NG matrix components. The expansion results in diagrams in (al), (bl-b4) and (cl-c2) 
in Fig (121) respectively. So we can see that our approach is very similar to the quark-quark scattering amplitude in 
Ref. Qjj if we open the condensate blobs. For example Fig. © (bl,b2) correspond to the quark-quark scattering 
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(c) (cl) 

FIG. 1: 2PI graphs up to two loops in NG basis. The solid and dashed lines are full propagators of quarks and gluons 
respectively. The black blobs at triple gluon vertices denote HDL. 



amplitudes with a one-loop Abelian correction to quark-gluon vertex attached to each quark line. Fig. ^ (b3) 
corresponds to the crossed box diagram for quark-quark scattering amplitude. Fig |(2])(b4) is proportional to <fi 3 and 
has no correspondence in Ref. [T^. One also sees that diagrams Fig. (cl,c2) are similar to Fig. @ (bl,b2) but 
contain one-loop triple-gluon corrections to the quark-gluon vertices. The main difference between our approach and 
Ref. 12| is that we work in super-phase while the authors of Ref. [12j worked in normal phase. One difference can 
be seen in the next section that there is a residue term in super-phase when deriving the Ward identity which has no 
counterpart in normal phase. Note that the difference of two approaches in normal and super phase has also been 
studied in Ref. 



IV. VERTEX CORRECTIONS 



In this section we will investigate the vertex contributions in the gap equation. An Abelian-likc Ward identity will 
also be derived explicitly from a cancellation between the Abelian and triple-gluon diagrams. 

There are two 1-loop diagrams which provide corrections to the quark-gluon vertex, one is the Abelian diagram 
and the other is from the triple gluon diagram, see Fig. ([3]). The full vertex is then the sum of two 1-loop diagrams 
and the tree-level vertex. Now we focus on the Abelian diagram Fig. (b) denoted by iA%. Note that the gluon 
line in zAi is a HDL-dressed propagator. Contracting zAi with momentum P we have 

iP„\\ a {P, K, K P) = g 3 J iD™ L (P')iT a ^iS(K - P')iT b P 

xiS(K - P - P')iT a j„ , (11) 

where iDf?„ DL (P') is HDL-dressed propagator. We note that f can be written as 

P = S-^K - P') -S-\K - P - P') . (12) 

Here we have neglected diagonal parts of selfenergy En and E 22 which is of order g<j> on the Fermi surface. The 
reason is that once En/22 ar e inserted into the vertex [the left full vertex in the first diagram of Fig. ([3])] in the 
gap equation whose right-hand-side is given by the first diagram of Fig. ([3]) , one will see that the contribution is at 
most of order <7 4 </>ln 2 <f> where the two logs are from the loop integral. We also assume that the gap is independent of 
momentum within the range \q — fj,\ < gfi around the Fermi surface, which means the gap is assumed to be constant 
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FIG. 2: 1PI graphs and their NG expansions. The solid and dashed lines are full propagators of quarks and gluons respectively. 
The triple gluon vertices are HDL-dressed. The blobs on quark lines are diquark condensates. 



in the gap equation with the momentum of the exchanged gluon being of order g\i. Inserting Eq. (p~2|) into Eq. (jlip . 
we see that S _1 (A - P - P') in Eq. (0 cancels S(A - P - P') in Eq. (TT])- But for S^A - P'), the procedure is 
a little more complicated. Because T b is not commutable with S(A" — P'), with a T b in the middle, S^-^A — P') in 
Eq. (fT2|) cannot directly touch S(A" — P') in Eq. (fTTj) . Making use of the commutator 



[T°, S _1 ] = - 



(T a J 3 
(T aT J 3 + J 3 T a )T 2l5 ^A e 







(13) 



we can move S _1 (A — P 1 ) over T b and cancels S(A — P'). Note that we have neglected the momentum dependence 
of the gap for the soft gluon exchange, therefore the above commutator is constant in the gap equation. Obviously 
the commutator only has off-diagonal elements which are of order <j>. After a short algebra, we obtain 



iP a K b l(T {P,K,K-P) 



g[iT,(K)T b - T 6 i£(A - P)} 
+gf abc T a {X nc (K) - S„ C (A - P)]T C 
d 4 P' 



D" DL (P')T a ^S(K - P')[T b , S" 1 ] 
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FIG. 3: Right hand side of the gap equation and the full vertex. 



(a) a al, jil, PI (b) ^ al, yd, PI 




K-P1-P2 K-P1-P3 
FIG. 4: HDL-resummed triple-gluon vertices 



xS{K - P - P')T a j v , (14) 

where S„ c is the quark selfenergy with the color part factorized out as X = T a S nc T a . 

Now we look at another one- loop diagram with triple-gluon vertex, Fig. © (c), which reads 

P a iA b 2a = g 3 P a J 0^iD™ L (PyD™ L (P'-P)iV v a %(P',-P,-(P'-P)) 

xT a i lu ,iS{K - P')T c i lt ,, . (15) 

The triple-gluon vertex iV^t is composed of two parts, the bare and HDL one iV = iV^ + iV HDL . First let us 
focus on the bare vertex as follows 

iy(°> = i V <P)F + lV (0)P (16) 

= -gf abc {[(-2P' + P) a9lw - 2P v9aii + 2P f ,g a „] 

+ [P'vQon + (P' - P)^9au}} , (17) 

where we have decomposed the bare vertex into a transverse and a longitudinal part. The first square bracket is the 
transverse part V^ F while the second one is the longitudinal part V^ p . Contracted with two HDL-propagators in 
Eq. (HU) V {0)p gives zero due to the fact that the HDL-propagator has the transverse property, while V^ F satisfies 
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a Ward identity as follows 



PriV^f = (P 2 -2P-P>)g^ 



(P' - P f - P 



9W 



(18) 



where we have factored out a constant —gf abc . Then we look at the HDL-resummed triple-gluon vertices as illustrated 
in Fig. (|3]). The first contribution is from Fig. (01(a) 



HDL 



(Pk 



Tr [ig 7fll T a HS {K x )i 91 ^ T a HS {K 2 )ig lll3 T° 3 iS (K)} 



^g 3 Tr[T ai T a2 T a2 ] J ^Tr [ 7w S (^i)7^S 



(K 2 ) lfl3 S (K)] 



(19) 



where K\ = K — Pi, K 2 = K — P\ — P 2 . The negative sign in the first line arises from the fermion loop. Note that we 
used the free quark propagator without condensate and sclfencrgy correction, because the condensate and selfenergy 
would contribute at higher order. There is a factor of 1/2 in the front due to the usage of NG basis. It is easier to 
work in normal basis and get rid of the factor 1/2. Hereafter we choose this way. Contracting iVj HDL with Pi Ml we 
obtain 



^i^™ aia2Q3 (l) = 5 3 Tr[T ai T a2 T a3 ] J ±^ 

x Tr [Pi So (#i ) 7m So (K 2 ) lll3 S (K)} 
We can rewrite Pi as Sq(K) — Sq (Ki) in the above formula and get 



(20) 



d 4 K 



x {Tr [S (K - P 1 ) 7fl2 S Q (K - P 1 - P 2 ) 7a13 ] 
- Tr [^ 2 S (K - P 1 - P 2 ) lm So(K)}} 
= igTr[T^T a ^} [U n ^ 3 (P 3 ) - H^ 3 (P 2 )] , (21) 

where n™^ 3 is the polarization tensor without the color part. Another diagram Fig. (J4Jd) is the same as Fig. ([4^) 
except that labels 2 and 3 are interchanged, which gives 



^i<M 2 M3 aiQ2Q3 ( 2 ) = fcnp^T^r*] [n- M3 (p 2 ) - rr^ 3 (P 3 )] 



Sum of two diagrams gives 

^^Cw 1 " 2 " 3 = ig^[T^{T a2 T^~T^T a2 )\ fc 3 (P 3 ) - (P 2 )] 

= -gf aia2a3 [n^ 2W (p 3 ) -n^ 2W (P 2 )] , (22) 

where n M2At3 is the polarization tensor with the color factor. We add Eq. (|22|) and (fT8|) together. By doing so, we 
make replacement a\a 2 a 3 — > bac, [i\[i 2 [i 3 — > avfi, Pi, P 2 , P3 — > P, —P 1 , P' — P, we obtain 

P iV abc — P \o\/(°) F '' abc _|_ jV HDL '> abc ] 

= -gf abc { [(P' - Pf - p' 2 ] 9»v - PWi" - P) - IV(P')]} 

= -9f abc [D HDL , v {P')-D H l DL ^{P'-P) . 
Substituting the above equation back into Eq. (| 1 5[) . we get 

P„iAl(K,K-P) = g 3 r bc J ^D™ L {P')D™ L {P' -P) 
x ~Dhdl„AP')- D hdl„A P '-P) 



(23) 



3 fabc 



t PI UP j_i 

xT a i lly nS{K - P')T c i 7 M 
d 4 P' 

(2^y 4 

xT a i7^iS(X-P')T c i7 ^ 



J0^[iD^(P'-P)-iD^(P')] 



(24) 
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Combining Eq. (fT4)) with ([24)) , we see that Eq. (|24|) cancels the second line in Eq. lfl"4)) . We obtain 



P a iA b a (K,K - P) = g[iH{K )T b - T b i£(if - P)} 

d 4 P' 



~9 



(2tt) 4 " v 



^DL ( p 0T a^ s(x _ p0[Tb;S -l ] 



xS(ir-P-P')T°7„ 



(25) 



Note that the last term results from the superphase because of the non-vanishing commutator [T & ,S _1 ] which is 
proportional to the diquark condensate. The similar term is there even in QED [43j |. We can add the bare vertex 
ig-f a T b to ik b a (K, K — P) and get the full vertex iY b a = igT b ^ a + ik b a which is the blob in the first diagram of Fig. 



PviT^ = ig[S-\K)T b -T b S-\K - P)] + I b x (K,K - P) 



(26) 



where I b x (K, K — P) with X = [S _1 , T b ] is just the last term in Eq. ([25)) . Here we have used Eq. Q and the property 
that Sq 1 is commutable with T 6 . 

With the identity in Eq. ((26)) we can evaluate the gauge dependent part or the £ part in the first diagram of Fig. [3) 

d 4 P 1 



<9< 



(2tt) 4 P 4 



[S^^T 6 -T 6 S _1 (^ - P) ~ igI h x{K,K - P)] S(K- P)T b lp P f 



p ■ 



(27) 



The first term inside square brackets is vanishing when sandwiched between on-shell wave functions. The second term 
is also zero due to / (PPPp/P 4 = 0. The contribut ion from l\ is evaluated in Appendix A and of sub-subleading 
order for gapped modes in the gap equation. 

Having P^A^ in Eq. (gSJ), can one derive ik b a ? In principle, one cannot. But at the limit P — > 0, one can derive 
the leading contribution of iA* . Note that there is a subtlety in defining the limit P — > because it is involved in 
two different types po — > 0, p — > and p — > 0,po — ► 0, which lead to different result for ik^ . We take the first limit 
for both sides of Eq. (|25|) to extract ik b as follows 



ik b = lim lim ik b AK,K~ P) ~ (r.h.s.) 
p^opo^o dpi 



(28) 



We consider the gapped modes. The contribution of I X (K, K — P) can be proved to be beyond the subleading order, 
see the Appendix A. Here we consider the first line in Eq. 



z.gT b ^--E(K - P) 
opt 







(29) 



P=Q 



which means both diagonal and off-diagonal parts are zero. We know that the diagonal parts S11/22 ~ <? 2 (&o — 
Po) ln(|fco — po\/M) which has no dependence on spatial momentum. The off-diagonal parts S12/21 is actually propor- 
tional to the diquark condensate which we assume in this paper do not have momentum dependence or the momentum 
dependence is of higher order. For gapped modes, the derivative of the second term of Eq. ([23)) will also give zero 
from the leading contribution, i.e. off-diagonal part which is of order g<p, see the Appendix A for the evaluation of 
the last term I X (K, K — P). We now take the second limit for both sides of Eq. |25)) to extract ik^: 



ik b 



lim \imik b (K,K - P) 

Po—>0 p— >0 



d 

dpo 



(r.h.s. 



The derivative of the first line in Eq. ([23)) gives 







igT" — £ (K - P) 
dp 



(30) 



(31) 



while the last term of Eq. ([25)) still gives zero. We can insert non-zero value of Ag into the gap equation where Ag 
couples to the Debye screened electric gluon whose contribution turns out to be of the sub-subleading order. 
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V. SUMMARY AND CONCLUSION 



We explicitly derived an Abelian-like Ward identity in color superconducting phase from Feynman diagrams, similar 
to the identity obtained by Nambu in normal superconductivity [43| . The identity arises from a cancellation between 
the Abelian diagram and the triple-gluon one. The same Ward identity was derived in Ref. [34[ in path integral 
approach. The identity has one additional term proportional to the Cooper condensate compared to that in normal 
phase, which is related to the gauge dependent part or £ part in one-loop vertex correction in color superconducting 
phase in Nambu-Gorkov basis. We finally reach that the vertex corrections are free of subleading contribution to the 
color superconducting gap by showing that the contribution of the additional term is beyond the subleading order for 
the gapped modes in the gap equation. The method proposed in this paper is equivalent to Ref. [l2j]. The difference 
between these two approaches is that our approach works in both super and normal phases while the approach of Ref. 
[l2l ] is based on the quark-quark scattering amplitude in the normal phase. One of the advantages of our approach 
is its compact form in NG basis, where a set of component diagrams including those with and without normal phase 
correspondences can be assembled into a single diagram in NG matrix form. 
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APPENDIX A: EVALUATE I b x (K, K - P) 
In this appendix we are going to evaluate l\{K, K — P) in Eq. ([26")) . We write it explicitly 

I b x (K,K-P) = - 3 3 J ^ DX U DL (PI)T^S(K - P>){T\S-i] 

xS{K - P - P')T: a lv . (Al) 

In 2SC, the commutator [T h , S _1 ] is given by Eq. (fl~3"|) . One only needs to estimate the order of 21- and 11-components 
of I x . The 11- and 22- components are of the same order. The 21-component is evaluated as 

(4) 2 i = J d*pi D ™ L (P') 

x(-T aT 7A1 ) {S 21 (i^)[T b , S- 1 ]i 3 S2i(K 2 ) + S 2 2(K 1 )[T b , S- x ] 21 S n {K 2 )} T a lv 



9 — [ d 4 P' D HDL (P') - - 

(2*)'J ar [ } kf-elkf-el 



= Citstb^A;^)^ / d 4 P'D^ DL (P') 



(2tt)< 

x(fc; + M ^fc 1 )(fc 2 °-M + fc 2 ) , n2 1 - — 2 1 - , (A2) 



jLO 2 _ ,2 U) 2 



1 c k!,i "-2 c k 2 j 

where we use notations K\ = K — P' and K 2 = K — P — P' and 

C X =C %XX = T aT J 3 (T b J 3 + J 3 T bT )J 3 T* = -^J 3 6 bs 

C 2 ,ij = T aT Li(T bT J 3 + J 3 T b )LjT a 

= nonzero, for b = 4, 5, 6, 7 and ij — 01, 10 
C 2 , 00 = for any b. (A3) 

We first consider C 2 ,n, i.e. gapped mode. We obtain 

= Ci^A^ / d^D^(P') 1 1 

\ 1 j /Cj e ki k 2 e k2 

x [{k° x +n- fci)(fc° - /i + hi) + (j) 2 } . (A4) 
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We take the derivative of the above with respect to modulus momentum p and then take the limit P — > as follows 

' [(k® + v - h)^" - n + k 2 ) + (f> 2 ] -2(fc 2 -A») 2k 2i k^ + fi-h 1 2k 2l 



(fcf-e*,) 2 fc 2 fcO 2 ^ fc o2_ £ 2 2 fc2 



x(-kii) 
= CiT2750(7 M A k 7 ; 
x(-kii) 



2(fci - M) 



(2tt)4 
(fci - /z) 



fc? 2 - ^ + 2^ 2 



d 4 F' D" DL (P') 



fei - A* 



fci - /j, 



2 \3 



where k\ = k\/ki. While the first term can be written as the derivative with respect to k: 

,) 2 



( d^P* D" DL {P>) ( 2 fcl M) = -1A [ d *P> D * DL (P>)^ 



i a 

29k 7 



a In 2 - + (3 In - + 7 
A* A* 



. 



Also the second term can be written as 

^-2^(4^-^)) 



/ d 4 P' D»J 



HDL (P') 



_9 
~5k 

d 



d 4p, D HDL {pl) 



2<t> 2 



(*? -O 2 



a In 2 - + /3 In - + 7 
A* A* 



Sk 4 
. 



2a In - +/3 
A* 



So that 



Use the same method, we can also get 



d(I b x h 



dp 1 
d(I b x hi 



= 



dp 



p=o 



p=o 



. 



p=o 



(A5) 



(A6) 



(A7) 
(A8) 

(A9) 



The contribution from C^io/oi (it is not vanishing) is zero for the gapped modes in the gap equation, which can 
be seen by 



C2,w/oiLiT a — 



(A10) 



So the only non- vanishing contribution from C 2l io/oi to the gap equation comes from the gapless modes, which is not 
bound. 

The 11-component of I b x is given by 



CVt.ToA^-J^ j d * P > D» V DL {P>) 

fcf - fi + ki 1 1 kl - n + k 2 

l,0 2 ,2 i,0 2 _ ,2 L.0 2 ,2 LO 2 _ ,2 

. K l e ki ^2 e k 2 ft l £ ki ft 2 e k 2 



(AH) 
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The gapless mode contribution is vanishing because 

T aT J 3 (T b J 3 + J 3 T bT )L T aT J 3 T b = 



It's easy to see that 



and 



T ai Lo ( T bi j s + J:iT b )J 3 T a J 3 T b = . (A12) 



P-0 dpi ^ J V '((fc?) 2 "^) 

- :L,AiL / D HDL (P') - ki 



2 Y d<f>J v ; (fc°) 2 -(e kl ) 5 

<7 2 k M , (A13) 



0(/£)ll/22 3 / jjHDL / pl\ 



((fc?) 2 -^) 2 ) 2 

~ 3 2 • (A14) 
Therefore we see that the contribution from vertices to the gap for the gapped modes is beyond subleading order. 
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